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SUMMARY 


1 A recurrence solution including the secondary effects of rotary inertia 
and shear deformation on discontinuous nonuniform beamlike structures is pre¬ 
sented for obtaining highly descriptive free-free natural mode data. Results 
of other studies are included to ascertain the significance of the secondary 
influences on the classical uniform beam. Numerical results are also included 
for the application of the method to a first-stage and fourth-stage configura¬ 
tion of a typical space research laionch vehicle. The n\unerical results indicate 
that shear deformation is generally the prime contributor of the secondary 
effects on typical launch vehicle configurations; the major influences of the 
secondary effects are seen in reductions in natural frequencies and in changes 
to the mode slopes, mode moments, and'mode shears. 


INTRODUCTION 


Natural vibration characteristics are frequently required for systems that 
can be appropriately represented as free-free beamlike structures. Such struc¬ 
tures are consistently encountered in the array of stages that comprise the 
typical launch vehicle for space research. 

An assessment of the natuiral vibration characteristics normally will 
include the calculations of the fundamental frequencies of oscillation along 
with a number of overtones. The deflected curves associated with each char¬ 
acteristic frequency, known as either mode shapes, characteristic functions, or 
eigenfunctions, are also desired data. In addition, the slopes, moments, and 
sometimes the shears associated with the mode shape are provided. Special 
definite integrals also are frequently computed, for example, the generalized 
or effective mass of a mode. 

These data have a variety of uses. Their most familiar use is in providing 
knowledge of natirral frequencies for designing in order to avoid states of 
resonant vibrations. This application is particularly important in designing 
or qualifying the spin program for unguided launch vehicles. (Trajectory dis¬ 
persion control is frequently accomplished by spinning a vehicle about its 
longitudinal axis.) Hi^ dynamic stresses and even structural failures might 


result if the roll frequency is near or coincident for an appreciable time to 
one of the natural frequencies of the system. 

In addition; the vibration response levels of a vehicle will generally be 
high near the frequencies of the natural vibrations of the basic structxire. 
Knowledge of the probable distribution of the disturbance spectrum can provide 
information for appropriate shock moimting of delicate payloads and can lead to 
a proper choice of the c3ynamic characteristics of instrumentation. 

Furthermore^ calculated mode shapes and frequencies can be useful in 
Interpreting measured vibration data. For example^ moment and shear data asso¬ 
ciated with calculated modal responses are frequently used for rapid assessment 
of the load significance of recorded flight data. 

Probably the greatest value of modal data is realized in series solutions 
of the differential equations of motion of the system for which they have been 
generated. The well-known orthogonality properties of mode shapes^ the fact 
that the modes satisfy the appropriate boundary conditions for their particular 
constraints, and the fact that the usual response of a beam system is adequately 
described by superposition of a few modes make them ideal functions for such 
applications, widely known as modal form solutions. 

The modal form series solution is particularly adaptable to formulating 
techniques for yielding gust, wind loads, ignition and separation responses, 
and other transient loads. The modal form approach has also proven popular for 
developing the characteristic equation of structures coupled with closed-loop 
autopilot systems. Such analyses permit investigations of the stability bound¬ 
aries of such systems and lead to proper gain levels and filter characteristics 
to yield stable performances. 

Such valuable applications of modal data have stimulated analysts to 
develop a variety of analytical techniques to achieve increased accuracy and 
scope of output with a minimization of input effort. Well-known and fundamental 
considerations of beam vibration are presented by Den Hartog in reference 1. 

The classical techniques of Rayleigh and Stodola are two of the methods dis¬ 
cussed in detail by Hartog. In reference 2, Scanlan and Rosenbaum describe the 
classical procedure attributed to Rayleigh and Ritz along with the Myklestad 
method. They proceed to meticulously outline a step-by-step matrix method 
employing influence coefficients. This latter procedure has been formulated 
into a rigorous matrix solution and is presented by Alley and Gerringer in ref¬ 
erence 3* Houbolt and Anderson have organized the method of Stodola and pre¬ 
sented other useful related material in reference 4 . An Interesting integral 
series solution technique that is also an effective solution to the nonuniform 
beam vibration problem is presented by Spector in reference 5- These references 
are but a few of the publicized methods relating to natural vibrations of non- 
uniform beams. 

When the analyst is confronted with a particular assignment for computing 
modal data, the choice of the most suitable method should be made in view of the 
objectives of the assignment and the virtues and limitations that characterize 
the methods. Energy methods following the concepts of Rayleigh or the ' 
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Rayleigh-Ritz approach generally yield adequate results on frequencies but 
possess weaknesses in the exactness and descriptiveness of the mode shapes and 
their derivatives. The classical Myklestad method as well as the matrix method 
of reference 5 treat the system as a number of liimped point masses coupled with 
massless springs. These discrete mass models also yield adequate frequency 
data but suffer a loss in descriptiveness of the other modal data, particularly 
in the shear ciirves associated with the modes. The Stodola process, or one 
similar, requires iteration to converge on both frequency and mode shape. 

Whereas the frequency is readily obtained to a desired accuracy, the mode shape 
is not so readily made everywhere convergent to a prescribed accuracy. In addi¬ 
tion, the procedure requires sweeping techniques to obtain modes of higher fre¬ 
quencies than the fundamental. This operation frequently presents practical 
problems due to degeneration in numerical accuracy. The integral series method, 
although capable of yielding both the fundamental and higher modes, presents a 
formidable computing problem, the maintenance of sufficient accuracy in per¬ 
forming numerically the required multi-integrations. In the original or classi¬ 
cal form, all of the aforementioned methods omit the secondary influences of 
rotary inertia and shear deformation. Many extensions to the classical forms 
have been made throughout industry and government, yet few satisfactory methods 
incorporating the secondary influences have been adequately described in gen¬ 
erally available literature. Shear deformation introduces an additional source 
of deflection to the customary flexural deformation considered in elementary 
beam theory. Cross-section rotary inertia provides additional dynamic loading 
to the system due to the rotational acceleration of the cross section of the 
beam. In most slender beams the shear and rotary inertia contributions to 
loading and deflection are small in comparison with those resulting directly 
from bending. In multistage launch vehicles of conventional fabrication the 
shear deformation is secondary to bending and the stioicture is generally so 
slender that rotary inertia may also be ignored. 

However, the use of fiber-glass stages has given unexpected emphasis to 
the importance of shear. The low ratio of shear modulus to the flexural modulus 
of elasticity of fiber glass increases the relative contribution of shear defor¬ 
mation to bending and thus produces significant shear effects in geometries 
that otherwise could be Investigated adequately without considering secondary 
effects. 

In addition, when computing the modal data of beam systems with low length- 
to-diameter ratios such as characterized by the upper stages of conventional 
multistage launch vehicles, the inclusion of both shear deformation and rotary 
inertia has been found to be important to data accuracies. 

This paper presents another technique for computing modal data on nonuni¬ 
form beams. The method has proven very satisfactory for accurately describing 
not only the mode shape, but its slope, mode moments, and mode shears. The 
method is inherently applicable to structures exhibiting numerous discontinu¬ 
ities in their mass and stiffness properties and requires no discretlng to an 
analogous lumped mass system. The differential equations of motion of the beam 
system are dealt with directly. Modes need not be computed consecutively as is 
required by a number of other methods. The technique has been employed suc¬ 
cessfully in obtaining accurate modal data up to the tenth mode on complex 
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structures. Furthermore, the formulation accurately incorporates the secondary 
influences of rotary Inertia and shear deformation. 

The organization of the problem in the matrix form that follows has proven 
to be hi^ly practical in obtaining numerical solutions on a high-speed digital 
computer. Also, the numerical integration technique that is an integral part 
of the method has been found to be hl^ly stable and accurate. 

The detailed mathematical developments of the method are presented with 
information pertinent to practical numerical solutions utilizing the digital 
computer. Some data are supplied for assessing the probable iD5)ortance of the 
secondary influences. An example of the application of the method to a typical 
four-stage launch vehicle is submitted with a study of the importance of sec¬ 
ondary influences on the upper stage. 


SYMBOLS 

p 

A cross-sectional area in shear, in. 


"A" 

solution 

for 

superposition method 




Al 

value of 


at xi 

due to a unit value of 

if 

at 

Xa, in.""^ 

A2 

value of 

V 

at 

due to a unit value of 


at 

Xa, Ib/rad 

a 

coefficient (eqs. (50)).> ^+1^^ Ih-rad^/in 




"B" 

solution 

for 

superposition method 




Bl 

value of 


at x^ 

due to a unit value of 


at 

Xa, rad/in.^ 

B2 

value of 

V 

at x^ 

due to a unit value of 

c 

at 

Xa, Ib/in. 

b 

coefficient ( 

[eqs. (30)), ^n+l^-> Ib-rad^/in, 

.2 



Cl 

value of 


at x^ 

due to a unit value of 

'Ir 

at 

Xg^, in./rad 

C2 

value of 


at x^ 

due to a unit value of 

if 

at 

Xa^ unitiess 

D 

maximum diameter of beam, in.; also used as determinant (eq. (52)) 

Bl 

value of 

C 

at x^ 

due to a unit value of 


at 

Xa^ unitless 

D2 

value of 


at x^ 

due to a unit value of 


at 

Xg, rad/in. 


k 


d 

E 

e 

f(a)) 

G 

I 

K 

k 

L 

M 

m 

nip 


^n+l 

P 

t 

V 

X 

X 

% 

Ax 

Z 


shear deformation coefficient (eqs. (jO)), Ih"^ 

modulus of elasticity in tending^ Ib/in.^ 

flexibility coefficient (eq.s. (jO))^ lb"^-in.“2 

frequency equation (eq. (37)) 
modulus of shearj Ib/in.^ 

area moment of inertia for a given cross section^ in,^ 
coefficient for cross section in shear 

effective spring constant of the rth mode; Ih/in, when mode 
shapes are considered dimensionless 

half increment (eqs, (30))^ Ax/2^ in. 

overall length of heam^ in. 

mode moment due to bending^ in-lb 

mass per inch of length, Ib-sec^/in.^ 

effective mass of the rth mode; Ib-sec^/in. when mode shapes 
are considered dimensionless 

matrix in recurrence solution of equation (27) 

point in derivation 
time, sec 

mode shear due to bending, Ib/in. 
column matrix of variables (see eq. (lO)) 
coordinate along length of vehicle, in. 

value of X at extreme left end of beam (lower boundary), in. 
value of X at extreme right end of beam (upper boundary), in. 
Increment in recurrence solution, " ^n^ 

rotary inertia, lb-sec^ 
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P matrix (see eq, (lO)) 

6 j_ matrix (see eq. (l 9 )) 

62 matrix (see eq. ( 20 )) 

^ deflection^ in. 

A matrix (see eq. (21)) 

A approximation of second derivative of variable (see eqs. (l2) 

to ( 15 ) and ( 21 )) 

f cross-sectional rotation due to bending^ rad 

00 circular frequency^ rad/sec 

increment of o) in trial solutions 

Subscripts: 

A "A" solution 

a lower boundary 

B "B’* solution 

I upper boundary 

n body station 

r mode 

Matrix notations: 

column matrix 

square dr rectangular matrix 

inverse matrix 
unit matrix 

Primed symbols denote differentiation with respect to x. 
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DERIVATION OF THE METHOD 


A derivation is presented for the differential equation of an oscillating 
beam structure expressed as a set of four first-order differential equations. 

A recurrence formula is established which relates the system variables (mode 
shapes, mode slopes, mode moments, and shears, and their derivatives) at one 
station of the beam to those at an adjoining station. Boundary equations are 
established for the free-free system and a characteristic equation is obtained 
for computing the critical frequencies of the system. Treatment of discontinu¬ 
ities is considered in detail, and the orthogonality and equilibrium relations 
are established. 


Equations of Motion 

A sketch is presented which shows the equations of motion. 
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Consider the differential element of a heam of dx length in the sketch, 
then the summation of the vertical forces on the element will yield 

„ ^ ax J dx . 0 (1) 

Summing moments about p and dropping second-order terms gives 

V dx + Z dx - ^ dx = 0 (2) 

From elementary heam theory, it is known that 


V = 



(3) 


M = El 


Sx 


(M 


where is the cross-section rotation due to bending which differs from 
because of shear deformation. 

Also, if interest is confined to the undamped rth natural mode, solutions 
of 'ilr, V, and M are assumed to have the form 


^(x,t) = ^r(x)e^‘^^ 

t(x,t) = i|ri.(x)e’'^^ 

V(x,t) = Vr(x)e^^^ 

M(x,t) = Mr(x)e^^^ 

J 

where is the angular frequency of londamped harmonic motion in the rth 

mode. 

It should also he noted that KA.G and El are functions of x only; 
therefore 

a(KAG) ^ d(KAG) 
bx dx 
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5(BI) ^ d(El) 

Sx ^ 

Substituting equations (5) in equations (l) to ( 4 -) and performing the indi¬ 
cated differentiations gives the basic form of the differential equations for 
the natural vibration of a beam, consideration being given to flexure, shear 
deformation, and rotary inertia. Henceforth the subscripts will be dropped for 
purposes of clarity, and the equations will appear as 


-mo^^ + — = 0 
dx 

(6) 

V - Za:^t - = 0 

dx 

( 7 ) 

^ +JL. t = o 

(8) 

dx KAG 


M dt _ 0 

( 9 ) 

El dx 


Expressed in matrix notation, equations (6) to ( 9 ) appear as 


0 -mo:^ 

0 0 

1 0 

0 0 


0 

0 

0 


0 10 
-Za? 0 1 


-1 



-10 0 0 


0 

0 


V 

1 

0 


c 

r 

0 

0 


V‘ 

0 -i- 

EI 

■. 

V 

M’ 

[mJ 


( 10 ) 


where the primes denote differentiation with respect to x. 

If the 4x8 matrix in equation (lO) is denoted as 3, and the 8x1 
column matrix of unknowns as X, then equation (lO) may be written as 


[6]{x> - 0 

The matrix equation (eq. (ll)) constitutes the equations of motion of the 
system. 
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Recurrence Solution 


The following relationships are introduced and are based on the assiuirption 
that the first derivatives of the system variables are linear over the small 
increment to (^n+1 “ ^n) 


dx dx 


+ AxAn(^) 




dx dx 


+ AxA^(ilf) 


W = ^ ^ An(t) 


dVn+1 ^ ^ 

dx dx 


+ AxAn(V) 


Vn+1 = ^ Ax + Vn + ^ An(V) 
dx 2 


dMp+i ^ 
dx dx 


+ AxAn(M) 


Mn +1 = ^ Ax + Mn + ^ An(M) 


At this point cetera^ are to be regarded as undetermined. 
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Equations (12) to ( 15 ) may be written in matrix notation: 



11 








The system variables must now satisfy the equations of motion at any specified 
station. That is, when x = x^+i. 
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by virtue of the constraint Imposed through the following determination of the 
A matrix. Substituting equation (22) in equation ( 23 ) gives 

3n+l 

Solving equation (24) for An yields 




r -] 


-1 

r- — 

[Pn+lJ 

62 


Pn+l 



(25) 


Thus, the variables A(0^ A(t), et cetera which were undetermined in equa¬ 

tions ( 12 ) to ( 15 ) are defined by requiring satisfaction of equation ( 23 ). 
Substituting equation ( 25 ) into equation (22) gives the recurrence formula for 
the functions V, and M: 




1 

H 

- 

62 


[pn+ij 

L. -J 

62 

-1 

Pn+l 


5 l 


L 












(26) 


It is expedient to perform the matrix manipulations to express equa¬ 
tion ( 26 ) in the form 



l4 


From equations (lO), (ll), and (20) it can be seen that 


HH 




For convenience, let 


0 


0 


0 




2 


-1 


a = mn+io:? 

t = 2.n+l<^ 


d = KAG 


-1 


n +1 > 


e = El 


-1 

n +1 


Ax 


2 

1 


KAGn+i 2 


0 


0 


1 Z^lX 


Eln+i 2 


(29) 


(50) 


Then^ equation ( 29 ) may be expressed as 



The determinant D of this matrix is given by the formula 


By utilising equations (lo), (ll), (19), (»), and (JO) and by perforMng 
the uSrlx operations Indicated In equation (26), the elements of the matrix 
P„tirSch Appears in equation ( 27 ), W be reduced to the explicit algebraic 

expressions given in eciuation (53). The elements are functions only of the 
sS^cSal properties at x^ and x^+i. of the interval length Ax, and of the 

unknown-freciuency o). Equations (32) and (55) provide the vorMng relationships 
for practical use of the recurrence formula given hy equation (2^). 


ik^[ek^(l-'bd)-^ 

ak[ek®(l-bd)-^ 

k 

1 

[ek2(l--bd)-^k 

[^k2(l-bd)-3 

ek2 

ek 

k(l+bek2) 

l+bek^ 

k2 

k 

[ek^(l-bd)-3k2 

[ek2(l-bd)-(3k 

ek5 

ek" 

aek^ 

aek^ 

ek^[ak^( l-bd)-'^ 

ek[ak®(l-bd)-'3 

ek2 

ek 

ek(l+adk2) 

e(l+adk^) 

aek** 

aek^ 

k(l+adk^) 

1+adk® 

ek5 

ekS 

ek^(l+adk2) 

ek(l+adk^) 

a(l+hek^)k 

a[1+bek®) 

ak2 

ak 

ak^ jek^(l'*bd)-^ 

ak jek^( 1-bd) 

aek5 

aek® 

a(l+bek2)k^ 

a(l+bek^)k 

ak5 

ak2 

k(l+bek2) 

l+bek® 

aek** 

aek^ 

ak^ 

ak 

(al^(l-bd)-'3 k 

[ak^d-'bd)-^ 

k 

1 

ek® {ak2{ l-bd) 

ek[ak®(l“bd) 

ak5 

ak^ 

[^k2(l-bd)-^k2 

[ak^(l-bd)-^ k 

k2 

k 

k^l+adk®) 

1+adk® 


(53) 


Boimdaiy Conciitions 

In the work that follows^ the solution will he restricted to considera¬ 
tions of the free-free vibration behavior. This mode of behavior is appropriate 
to the in-flight characteristics of launch vehicles. Henceforth, the sub¬ 
scripts a and Z indicate that the associated quantity is evaluated at the 
lower boundary x = Xg^ or at the upper boundary x = x^^ respectively. The 

boundary conditions at both ends of a free-free beam require that 
Condition (a): 

Va = = 0 

Condition (b): 


Ma = = 0 
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Substituting condition (a) in equation (8) shows that 
Condition (c); 


Condition (d) 


^ 




Substituting condition (a) in equation (7) gives 


Condition (e): 


Condition (f): 


K = -2-a^ta 


= -Zia^n 


Substituting condition (b) in equation ( 9 ) yields 


Condition (g) 


ia = tz = 0 


Solving equation (6) for dV/dx gives 


Condition (h) 


Condition (i): 




v; . mjopt 


These conditions may now he expressed, in matrix form. For the lower 
boundary conditions it follows that 
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r 






^a 


1 

0 



Ca 


0 

1 

< 


0 


0 

0 



ta 

1 = 

1 

0 

< 

“a'^^a 


0 


Va 


0 


0 

0 

K 


-Za^ta 


-Za^ 

0 

Ma 


0 


0 

0 

> 


^ J 






and for the upper boundary conditions 


r 








h 


1 

0 

Cz 


Cz 


0 

1 

'1^1 


0 


0 

0 

n 

>- < 

tz 

i 

1 

0 

VI 


= 

0 

nizO^ 

Vz 


0 


0 

0 

Ml 




-Zz^ 

0 

Mz 


0 


0 

0 











( 5 ^) 


(35) 


Calculation of Frequencies 

Suppose now that a number of stations have been established at definite 
points along the span of the beam. The system variables at all stations 

are then uniquely determined if values are assigned to 
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Equation {^k) determines the system variables at x = Xa- Equation (27) yields 
the system variables at x = Xj^ in terms of the values at x = Xa« Once the 
variables at x = Xg^ are knovm, equation ( 27 ) may be used again to compute the 
variables at x = and this process may be continued until the system 

variables at all stations are known. This process will be termed "a solution 
of the system" for the given values of tg, and co. The principle of 

superposition of solutions is employed by using an arbitrary value of 01 in 
two separate solutions. The first solution has ^g = 1, ta = 0 and is called 
the "A" solutionj the second has = 0, 'ifa = 1 and is called the "B" solu¬ 
tion. Then by using conditions (a) and (g) of the section on bo\indary condi¬ 
tions, it can be concluded that 



The subscripts A and B designate the separate solutions for the afore¬ 
mentioned boundary constraints. The coefficients A^^, A 2 , Bj_, and B 2 are 

essentially influence coefficients where A^^ and are the values of 

and V, respectively, at x^ due to a unit value of if at Xg. The values 
of if* and V, respectively, at x^ due to a unit value of ^ at x^ are 

termed B 2 _ arid B 2 . 

Nontrivial solutions of equation ( 56 ) exist if, and only if, the determi¬ 
nant of the square matrix is zero. The natural frequencies are therefore deter¬ 
mined by the equation 


Al 

A2 



(57) 


hereinafter referred to as the frequency equation f(co). The coefficients of 
equation ( 37 ) are frequency dependent as is evident from inspection of the 
Pjj+2. matrix defined by equations (jO), (52), and ( 55 ). This frequency depend¬ 
ency permits the satisfaction of equation ( 57 ) and yields the natural frequen¬ 
cies of the system. Also, from superposition of the "A" and "B" solutions, the 
following additional relationships can be stated: 
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(38) 



where C]_, C 2 , and D 2 are influence coefficients for ^ 7 * 

From the second row of equation ( 36 ) 


Hence, 


A2ta + B2^a = 0 



From the first row of equation ( 58 ) 


(39) 


(40) 


= 


Solving equation (4l) for 

gives 


Cll^a + I^l^a 

and substituting equation (4o) for 


(41) 


'I'a 


■^2^1 " B2C 


and normalizing in tenns of = 1 yields 


t = ^2 

A2D1 - B2C1 


Substituting equation (42) in equation (40) yields 


(42) 


ta 


-B2 

A2D2_ — 22^2^ 


(43) 


It is reminded that the appropriate values of A;]_^ k 2 ) B 2 ^ 

C 2 , Di, and D 2 are only those associated with "A" and "B" solutions 
satisfying equation (37)" The numerical procedtire in obtaining valid solutions 
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is to perform trial solutions for assumed values of cjd until equation ( 37 ) is 
satisfied. The lowest frequency satisfying the equation yields the fundamental 
natural mode of vibration. Succeeding higher frequencies, which satisfy equa¬ 
tion ( 57 )/ yield the progression of overtones or higher modes. 


Mode Shapes and Related Data 

The coefficients to D 2 associated with the proper values of cd 

which satisfy equation (37)^ when substituted into equations (i^2) and (^3), 
yield appropriate initial values of ^a mode shapes. The 

recurrence formula of equation ( 27 ) relates the modal characteristics at a given 
station to an adjacent station. Thus, having the initial values leads to a 
complete solution for all stations. Such a solution yields the mode shapes, 
cross-section rotations, mode moments, mode shears, and the first derivatives 
with respect to x of each function. 

It should be noted that in letting t^i equal imity in obtaining equa¬ 
tions (4-2) and (^<-3), the system variables derived from use of the two relation¬ 
ships are henceforth compatible with the mode shape normalized at x = x^ 

(i.e., ^ 2 = 1 ). For purposes of preventing conflicts in units, it should be 

understood, that the form of the system variables, when = 1, is in essence 



m 


That is^ V, and M from a final solution are actually the ratios 

respectively. These ratios are the system 

variables per unit of mode deflection at x = . 


Treatment of Discontinuities 

In most practical problems associated with nonuniform beams, El, KAG, 
Z, and m exhibit many discontinuities over the range Xg^ ^ x ^ x^ . Conse¬ 
quently, provisions must be made for a convenient solution across discontinu¬ 
ities. When a discontinuity is crossed, the system variables have to be 
reevaluated. 

From physical consideration >(f, M, and Y are continuous, and it 

will be shown that M*, and V’ are discontinuous. 

Let the subscripts (-) and (+) denote values of system variables at a 
discontinuity which is approached throu^ values of x lower than or higher 
than, respectively, the value of x at the discontinuity. 
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Then from equation ( 9 ) where M(_) = M(+) 




dx 


EI( + ) 


From equation (8) where V(_) ~ ^( + ) ~ '^(+) 


(^5) 


and from equation ( 7 ) 


4t(.) 


dx 


= i(.) 


(-) 


KAG 


( + ) 


dM 

“dx 




(116) 


(^7) 


From equation (6) with “ ^(+) 


dV 


(+) 


dx 


= m, 


(+)'^ 9 -) 


(>*8) 


The relationship across a discontinuity may then he shown in matrix form 



0 + 


0 

0 

1 

0 

0 

'( + )' 
0 


•Z/ . \CjD? 


KAG(+) 

0 

0 



(i^9) 
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This relationship is exactly satisfied when /Xx: - 0 is entered into Pn+l 

of equation (55). This contributes to the efficiency of the method for use in 
conjunction with high-speed digital computers because the same recurrence for¬ 
mulas are used to relate system variables across a discontinuity as are used to 
relate variables across a finite interval. 


Orthogonality Condition 

One of the major uses of natural mode data is in modal form series solu¬ 
tions. Their practical applications in series solutions is a result of the 
many simplifications that result from their orthogonality relationships. These 
relationships for modal solutions, restricted to elementary beam theory, are 
well known and widely referenced in reports and standard texts on elementary 
vibrations. The addition of rotary Inertia and shear deformation alters the 
conventional relationship and a clear description of the orthogonality rela¬ 
tionship for nonuniform free-free beams with secondary effects is considered 
essential to the completeness of this paper. In reference 6 , Leonard derives 
the orthogonality relationships for the uniform free-free beam with shear 
deformation and rotary inertia. In the following discussion, the conditions 
for the nonuniform free-free beam will be set forth. 

Substituting equation ( 8 ) into equation ( 6 ) and adding subscripts to ^ 
and ^ 1 ^ to designate the specific mode r yields 


. CL 

““i- to 


/dU ' 


= 0 


(50) 


Also, substituting equations ( 8 ) and (9) in equation (7) gives 




= 0 


(51) 


Multiplying equation (50) ( 5 l) adding both results 

together, and integrating over the length of the heam, yields 




Xa 


to 




to 


X 


'l^n El — Idx 


X 


a 


to\ dx 

d^r \ 

KAGlet'll ^ 


(52) 
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since the r and n subscripts have been selected arbitrarily, equa¬ 
tion (52) is equally valid when expressed by reversing r and n. Making 
this reversal yields 



Expanding equations (52) and (55) by integrating, by parts, the first and 
second terms on the rl^t-hand side of the equations and then subtracting the 
extended version of (55) from the extended form of (52) will yield 




dx 






X- 


- tnEI 


dt] 

dx 


<iin 

+ tpEI — 
Y dx 

^a 


^a 


( 5 ^^) 


For the free-free boundary conditions 


KAG^^ - tn 




X- 


= Vn(xa) - Vn(xj) = 0 


El 


d^ 

dx 




X, 


K Equally valid 
' for n or r 


= Mj^(x^) - Mn(xa) = 0 


(55) 


y 


With these relationships^ equation ( 5 ^) leads to the final conclusion that 



- 


)/""^(“^n^r + ZVr)'3^ = 0 


(56) 


2 k 


and when ^ it is necessary that 



(mC„tr + Z*n*r)'^ = « 


(% / %• (57) 


When (% = 0^^ it is usual to define the preceding integral as the effec¬ 
tive mass of the mode. 


mr 




(% = 0^ (58) 


To avoid confusion and to obtain the effective mass in tenns of mass units 
(lb-sec2/in. it is necessary only to interpret the symbols and as 

applicable to the dimensionless mode. 

Now^ instead of subtracting as was done to obtain eq^uation (5 ^)} "tbe equa¬ 
tions are added; then with the additional results of equation (5?) when n ^ r, 
it is found that 



+ kag(^ 

dx dx \dx 



n / r (59) 


and when n = r, it is usual to define the preceding integral as the effective 
spring constant of the mode: 



where the same statement regarding units for equation (58) applies. 

Equations ( 57 ) ( 59 ) s-re the orthogonality properties of the free-free 

beam of variable mass and stiffness with shear deformation and rotary inertia. 

Equations (58) and (60) are useful relationships defining the effective 
mass and effective spring of the rth mode referenced at the normalizing point 
of the mode shape ^3- = 1. 
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Equilibrium Relationships 


In addition to the useful orthogonality relationships, two other valuable 
integrals associated with the eq\iillhriuim of the mode are frequently encoun¬ 
tered. For the free-free system to possess dynamic equilibrium while under¬ 
going natural vibrations, the following relationships prevail: 



m^P dx = 0 


(61) 




( 62 ) 


SOME ASPECTS OF APPLICATIONS 


Basic relationships for the recurrence solution of the beam vibration 
problem are given in equations (27)^ (jO) ^ (32); and ( 33 )- "^tie natural fre¬ 
quencies are obtainable from equation ( 37 )« Orthogonality and equilibrium 
relationships are given in equations (57)^ (59)^ (61); and (62). 


Descriptiveness of Output Data 

The type of solution set forth by the foregoing equations is inherently 
adaptable to producing highly descriptive modal data. It has been shown by 
equation (49) that ^ and Y' are discontinuous functions by vir¬ 

tue of discontinuities in El, KAG, Z, and m. The quantities obtained by 
the solution are of equal descriptiveness to that of the input data; that is, 
for every discontinuity in input data (EI, KAG, Z, and m) there will result 
an associated discontinuity in the output functions. The quality of these data 
cannot be equaled by the discrete mass methods of analysis such as that of ref¬ 
erence 3- The method is also inherently suitable for automatic digital¬ 
plotting routines since the input can be supplied or generated by interpolation 
to a definition as fine as desired to give an acceptable plot of output. 


Numerical Solutions 

The successful use of the foregoing formulation is largely dependent on 
the details of the computer program. Several important features to be observed 
are set forth in the following text. 

Machine time .- Improper programing can result in lengthy computing times 
yielding unreasonably costly data. Conversely, proper programing for minimum 
computer time can be most rewarding with the subject formulation. For example. 
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programing the inverted and extended form of equation (28) given by equa¬ 
tions (32) and (33) reduced the machine time to one-fourth that required by a 
program performing internally the operations of equation (28). Unduly stringent 
accuracy requirements in determining the natural frequencies by equation (37) 
can also result in exorbitant machine time. Extensive use of the method of 
this paper has indicated that three modes with frequency convergence to 0.1 per¬ 
cent error can be achieved with around 3OO stations of input data to an IBM 7 ^ 9 ^ 
electronic data processing system in about 5 minutes. It should be noted that 
this accuracy on frequency is valid for the mathematical model but it is not 
indicative of the accuracy of predicting the frequencies of the actual struc¬ 
ture. The accuracy is obviously dependent on the appropriateness of the math 
model. 

Attempting to satisfy upper boundary conditions to unneccessary accuracies 
can also be costly. For exan^le^ and are zero for the free-free 

boundary case. These absolute conditions cannot be achieved by digital com¬ 
puters. The analyst^ however^ should accept finite boundaiy values that are a 
fraction of 1 percent of the peak absolute values of their respective functions 
over the total span ^ x ^ X2. Boundary values which are zero correct to 

the fifth significant figure of the maximum value of a f;inction can readily be 
obtained. 

Eear zero stiffness .- In applying the recurrence technique of this report 
the unwary analyst can be confused and misled by using inputs of KAG and El 
that are near zero at or near the upper boundary position. This condition is 
encountered particularly with pointed nose cones or other similar structures 
that essentially taper to zero. The problem manifests itself in a radical 
variation of the mode shape in the area approaching the upper boundary of the 
span. This variation is a combined result of near zero stiffness of KAG 
and/or El and the failure to achieve numerically the absolute theoretical 
zero boundary condition. This condition is readily corrected by avoiding^ near 
the free ends, KAG and El input values of less than 0.0001 of their respec¬ 
tive average values over the total span. 

Superposition for final modes .- The sequential nature of a solution by 
the method proposed in this paper is as follows: Firsts obtain for a variety 
of frequencies the influence coefficients Aq^ A2^ and B2 by the "A" 

and solutions discussed in "Calculation of Frequencies." From this trial 
technique^ a frequency is discerned to a desirable accuracy that will satisfy 
the frequency equation (eq. ( 37 ))* Second, after obtaining a natural frequency, 
the influence coefficients associated with the critical frequency are substi¬ 
tuted into equations ( 42 ) and ( 43 ) to obtain the boundary values and ta 

associated with the normalized natural mode. Finally, with the knowledge of 
0), and initial conditions and coefficients can be fully ascer¬ 

tained, and the recurrence solution of equation (27) can be extended over the 
total span. This solution yields data on all of the system variables for all 
desired stations. An important feature associated with this last operation 
should be noted. Entering and ^Ifg simultaneously into equation (27) in 

the final solution for the modal characteristics is a different digital opera¬ 
tion than the original procedure of superposition that led to knowledge of ^g 
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and Consequently^ a simultaneous solution will result in different dig¬ 

ital round-off errors that will reduce the accuracy of the final boundary 
values below that previously established by superposition. Since a large niun- 
ber of significant figures are essential to the mechanism of the numerical 
solution^ a duplicate procedure for the initial operation is recommended in 
obtaining the final modal data. The available values of and ta should 

be used in independent ‘’A" and solutions and should be combined at corre¬ 
sponding X stations to obtain the desired final modal data. 

Integration intervals .- The fundamental differential equations of motion 
given by equation (10) ar^ all of first order. Integration is performed on all 
variables with the assumption of linear variations of the first derivatives as 
defined by equations (l2) to ( 15 ). For this assumption to be acceptable, the 
Ax interval must at all stations be suitably small. Studies on interval sizes 
have shown that adequate results are normally obtainable by describing the 
physical characteristics of a vehicle at all points of discontinuity with the 
added constraint that Ax = L/IOO. The latter constraint is easily programed 
into the computer and obviates a large amo\int of repetitive input data over 
long constant sections. In addition, where continuous but nevertheless radical 
variations occur in the input functions, additional stations should be included 
to Insure that the variations are adequately defined. 

Frequency search technique .- The economy of the machine solution is inti¬ 
mately associated with the iterative technique employed to obtain the unique 
frequencies that will satisfy equation (37)• Equation (37) is of the form 

f(a)) = 0 

and must be solved by trial and error. 

Successful results have been achieved by performing repeated trial solu¬ 
tions for f(a)) by systematically increasing the frequency cd by Aco -until a 
sign change is indicated, i.e., 

t{Gp) 

where aP is the nth trial value. At this point, a second-order (parabolic) 
curve fit is applied to points at and and the 

resulting analytical expression can be solved for o) where f(a)) = 0. This 
operation can be continued with each improved frequency until the difference 
between successive frequencies is less than some preassigned tolerance. 

Convergence can be hurried by good estimates of for the specific 

problem and refraining from a too stringent tolerance in frequency. Increasing 
Ao) and the frequency tolerance for successive overtones is also advantageous 
in maintaining machine time within practical limits. 
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Secondary Influences 


A principal olDjective of the mathematical development of this paper was to 
provide a means for including^ as well as investigating^ the effects of the sec¬ 
ondary contributions of shear deformation and rotary inertia in modal analyses. 

Comparison with classical cases .- Studies have been made by comparing the 
classical or exact modal data for uniform free-free beams without secondary 
effects with data obtained by the subject recurrence technique with secondary 
effects. The two methods, secondary effects being excluded, agreed to an accu¬ 
racy of greater than 0.01 percent on the first three mode frequencies, mode 
shapes, and nodal points. It was felt that this agreement qualified the numer¬ 
ical technique for use in further comparative studies. 


The salient results of some studies to determine the significance of shear 
and rotary inertia on the classic uniform beam are furnished in figure 1 . Data 
were derived for a solid cylin¬ 
drical beam and for a thin-wall Length-to-dlameter ratio, l/d 


cylindrical beam of a thickness- 
to-dlameter ratio of O.OI 5 . 

Since the effectiveness of shear 
deformation is dependent upon 
the ratio of the shear modulus 
G to the bending modulus E, 
the studies were made for both 
fiber-glass and steel beams. It 
can be seen from the curves for 
the steel beam, that for either 
of the cross sections the com¬ 
bined secondary influences 
caused frequency reductions in 
the first three modes of not 
greater than 5 percent when 
ratios of length to diameter 
(L/D) are greater than 20 . This 
frequency reduction compares 
with those reductions of up to 
l 4 percent for values of L/D 
greater than 20 for fiber-glass 
beams. The value of G/E for 
the fiber-glass beam was 0.20 
whereas that for the steel beam 
was 0 . 42 . It is also evident 
that the frequency reductions 
for both the steel and fiber¬ 
glass beams are more pronoTmced 
for the thin-wall cases than for 
the solid beams. As the aspect 
ratio (l/d) decreases, the reduction in 



(a) Steel beam. 

Length-to-diameter ratio, l/d 



(b) Fiber-glass beam. 


Figure 1.- Effects of rotary inertia and 
shear deformation on the frequencies 
of free-free vibrations of uniform 
cylindrical beams. 

frequency becomes quite pronounced and 


reaches more than 50 percent in the third mode for the thin-wall cases for 
aspect ratios of 4 and 7 for the steel and fiber-glass beams, respectively. 



In reference 'J, Kruszewski presented closed-form solutions for the uni¬ 
form beam -with shear deformation and rotary inertia. The report gives general 
data for readily determining the significances of the secondary Influences for 
general beam cross section and elastic properties. The curves of figure 1 
obtained by the numerical procedure of this paper were checked against the data 
from the closed-form solution of reference 7 and were found to be in nearly 
perfect agreement. This coinparison further qualifies the numerical technique 
as a satisfactory and accurate procedure for incorporating the secondary influ¬ 
ences into the beam solution. 

Effects on launch vehicle .- The six data points given in the fiber-glass 
curves of figure 1 show the frequency reductions due to secondary influences on 
a four-stage space vehicle with L/D = 25 and the fourth stage of a space 
vehicle with L/D = 5 . Of the first-stage length 5.5 percent was fiber glass 
and 26 percent of the fourth-stage length was fiber glass. The first-stage data 
reveal reductions greater than those predicted for thin-wall steel beams but 
are less, as would be expected, than those for thin-wall fiber-glass beams. 

The fourth-stage frequencies for L/D = 5 were significantly reduced by 
the secondary effects and closely correspond to the reductions predicted for the 
thin-wall fiber-glass beam. It appears that the fiber glass appreciably con¬ 
tributed to the frequency reduction even though the fiber glass extended over 
only one-fourth of the stage length. The fiber-glass section extended over a 
span running from 0.25 to 0.55 of the length which was a region subjected to a 
significant proportion of the shearing action. The specific vehicle for which 
these data were derived is documented further in the section entitled 
"NUMERICAL EXAMPLE." 

Measure of significance .- The aforementioned considerations suggest the 
usefulness of the curves of figure 1 in estimating the probable sl^ificance of 
shear defomiation and rotary inertia in free-free natural frequencies of launch 
vehicles. Both the vehicle data and the curves clearly Indicate the necessity 
of including secondary influences if accurate analyses are desired on low-aspect 
ratio structures. 


NUMERICAL EXAMPLE 


An example is presented of an application of the recurrence solution to 
the actual research vehicle illustrated in figure 2 . Solutions have been pro¬ 
vided for the first- and 


Eirst-stage separation Second-stage separation Third-stage separation 



Figure 2.- Typical multistage research 
launch vehicle. 


fourth-stage configurations 
of flight. These two cases 
are submitted to show an 
actual numerical example 
derived from the subject 
method and to Illustrate 
the varying importance of 
secondary effects with the 
vehicle * s length-to-dlameter 
ratio. The basic physical 
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characteristics of the vehicle required as input to the program have been 
recorded and submitted in table I. The tabulation provides input for both 
cases by taking the physical characteristics at x = 619.85 as the first 
applicable quantities for the fourth-stage analysis. 


Frequency Data 

In table II, the frequencies and percent reduction in frequencies are 
given for the first three modes of vibration for the two stages of the vehicle 
of figure 2. The data are displayed for analyses without the secondary effects, 
with both effects, with rotary inertia only, and with shear deformation only. 

Inspection of the data for the first stage will show only small effects of 
shear deformation and rotary inertia. The influences grow progressively with 
increasing modes as would be normally anticipated. The maximum contribution of 
rotary inertia and shear deformation is seen to be a reduction in frequency of 
slightly less than 5 percent on the third mode. The first-stage configuration 
has a length-to-diameter ratio of approximately 25* In the fourth-stage con¬ 
figuration with an aspect ratio of 5.0, the importance of rotary inertia and 
shear deformation are quite evident. The combined secondary influences produce 
a reduction of 58 percent in values computed without consideration of shear 
deformation and rotary Inertia. The significance also increases with increasing 
modes. These results strongly suggest that for reliable frequency calcalations 
on the upper stages of typical launch vehicles, consideration of the secondary 
contributions to flexure must be given. 

It is interesting to note, however, that good results would have been 
obtained in all cases shown by incorporating in the solution the effects of 
shear deformation only. The frequency error most affected by the secondary 
influence (third mode of the fourth stage) woiild be slightly less than 7 Per¬ 
cent if shear deformation only had been considered. 


Elastic Curve Characteristics 

A graphical comparison of the modal functions for the first three natural 
modes for the first stage of the example vehicle would reveal only trivial 
departure from concurrency. However, it will be shown later that this is not 
the case for the fourth-stage configuration where the secondary effects have 
been shown to be appreciable. 

Mode shapes .- In figure 5; the fourth-stage mode shapes ^ are compared 
for solutions with and without secondary effects. The comparison of mode 
shapes was accomplished by the application of the method of least squares. In 
spite of the severe reductions noted in frequencies, the first three mode 
shapes retain similar graphical characteristics. 

Mode shape inflections near tips .- An interesting elastic curve character¬ 
istic should be noted that results from the inclusion of shear deformation. 

For beams with secondary effects, an inflection in the mode shape is frequently 
observed near one or both of the free ends and outboard of the outermost nodal 
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2 - Without rotary inertia and shear deformation 

[* -With rotary inertia and shear deformation 




X, inches 


Figure 3.- Effects of rotary inertia 
and shear deformation on mode 
shapes of a low-aspect-ratio 
upper stage of a launch vehicle. 


point. Inspection of either of the 
three inodes of figure 5 'wlH indicate 
the inflection or reversal in curva¬ 
ture on the right-hand end in the 
vicinity of x = 90 inches. The 
inflection point occurs vhen the sec¬ 
ond derivative of the mode shape 
hecoines zero. An expression is 
readily obtained by differentiating 
ecpiation ( 8 ) that indicates the 
parameters involved in the phenomenon, 
that is 


since 


d 2 ^ d / V \ ^ 

^ j " 



(65) 


at the Inflection point, and 



< 0 


to the ri^t of the Inflection point. Then it follows from equation (65) that 
for the dip in the elastic curve on the end to occur 


d__v_>^ (64) 

dx KAG dx 


Equation (64) simply states that when curvature contribution due to shear 
deformation is greater than the curvature contribution to the elastic curve due 
to bending, the inflection near the tip will be evident. In elementary beam 
solutions, shear deformation is not considered) this case corresponds to the 
limiting case where KA.G —> <». It is then obvious from equation (64) that for 
cases where shear deformation is ignored, the shear deformation contribution to 
curvature is zero and the conditions of equation (64) can never be realized. 


Mode slopes.- A pronotmced departure between data computed with and with¬ 
out secondary effects is seen in the comparison of mode slopes in figure 4. 

The solid curves are the slopes computed without shear deformation and rotary 
inertia. The dashed curves were computed with secondary effects. The slopes 
are compatible with the amplitudes of the mode shapes given in figure 5 
established by the least-squares method. The most conspicuous differences are 
the discontinuities that are evident in the data incorporating rotary Inertia 
arii^ shear deformation. The discontinuities result from shear deformation and 
are not evident in results without shear deformation, whether Including rotary 
inertia or not. Accurate mode slopes frequently are required in system 
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stability studies of autopilot 
structural feedback. It is evident 
from figure k that, for accurate 
structural feedback analyses on 
short L/D configurations, the sec¬ 
ondary effects, particularly shear 
deformation, should be included in 
developing the req^uired mode slopes. 


Mode slopes compared with 
cross-section rotation .- In elemen¬ 
tary beam theory, which ignores both 
rotary inertia and shear deforma¬ 
tion, the cross-section rotation is 
quantitatively the sajne as the 
slope. The differences between the 
slope of the elastic curve ^ * and 
the cross-section rotation ^ are a 
further measure of the significance 
of secondary effects. A comparison 
of these modal fiinctions is shown 
graphically in figure 5 ^or the 
fourth stage of the example vehicle. 
The observed differences between 
^ * and ^ are completely a 
result of shear deformation and 
are not evident in solutions 
ignoring shear. 


Load Characteristics 



X, inches 


Figxire 4.- Effects of rotary inertia 
and shear deformation on mode 
slopes of a low-aspect-ratio 
upper stage of a launch vehicle. 



of Modes 


The moment and shear dis- - 

tributions for the fourth stage ^ 

of the example vehicle of fig- X, Inches 

ure 2 are submitted in figures 6 

and 7* Comparative curves are Figure 5-- Comparison of cross-section 

given for data calculated with rotations and mode slopes of a low- 

and without the secondary influ- vehicle. (All Lrves include the 

ences of shear deformation and effects of shear and rotary 

rotary inertia. The moment and inertia.) 

shear values are compatible with 

the relative mode amplitudes shown 

in figure 5 s-s established 

by the least-squares method of comparison. Both the moment and shear data show 
significant differences with increasing modes between solutions with and with¬ 
out secondary influences, the magnitudes of the differences increasing with 
increasing modes. An appreciable separation in the right nodal point of the 
third mode moment curves is observed. The startling differences that are 
observed between amplitudes of the comparative data for the second and third 


modes are principally attributed to the lower frequencies and consequently 
lower inertia loading on the beam for the solutions with secondary effects. 




_1 _!- -J-i-' 

20 ho 6o 8o 100 

X, inches 


Figure 6.- Comparison of "the mode 
moments from solutions with 
and without secondary effects 
for a low-aspect-ratio upper 
stage of a launch vehicle. 


Figure Comparison of mode shears 
from solutions with and without 
secondary effects for a low- 
aspect-ratio upper stage of a 
launch vehicle. 


CONCLUDING REMARKS 


A recurrence solution is presented that is especially suitable for 
obtaining highly descriptive modal data on severely discontinuous nonuniform 
beamlike structuresj the solution includes the secondary influences of rotary 
Inertia and shear deformation. Data are provided for helping the analyst in 
estimating the probable significance of shear deformation and rotary inertia 
in free-free natural frequencies of launch vehicles. 

Numerical examples of the solution are included for a first-stage and a 
fourth-stage configuration of a research vehicle. The secondary influences of 
rotary inertia and shear deformation had trivial effects on the first-stage 
data, but were found to have very significant effects on the data for the 
fourth stage. 

In most practical applications to launch vehicles, it is concluded that 
shear deformation is considerably more important than rotary inertia. Accept¬ 
able results could generally be obtained by including shear deformation only. 
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The authors’ experiences indicate that^ in general^ for most beam applica¬ 
tions the secondary influences have the least obvious effect on the mode shapes. 
An appreciable influence is noted in reducing frequencies and large effects are 
experienced on the mode shears and moments. The appearance of discontinuities 
in the mode slopes is a further influence resulting from the consideration of 
shear deformation. 

The addition of shear deformation to elementary beam vibration theory can 
result in an inflection in the mode shape outboard of the outermost nodal 
points. 

Successful and economical solutions by the subject method are intimately 
dependent upon details of the computing program. Some salient points are dis¬ 
cussed on appropriateness of input data^ superposition problems^ numerical 
integration intervals^ and frequency iteration routine. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April 9; 1965* 
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TABLE I.- PHYSICAL CHARACTERISTICS OF A MULTISTAGE LAUNCH VEHICLE 


in. 

lb-sec2 

lb-sec2/in.2 

lb-in. 2 

KAG, 

lb 

-5'^.35 

1.70 

0.01743 

0.100 X io9 

10.00 X 10^ 

-19.35 

1.70 

.01743 

.100 

10.00 

-19-35 

14.70 

.06030 

14.770 

33.54 

0 

14.70 

.06030 

14.770 

33.54 

17.85 

14.70 

.06030 

6.745 

35.54 

17.85 

14.70 

.06030 

6.745 

53.54 

19.35 

14.70 

.06030 

6.745 

35.54 

19.35 

13.39 

.13094 

6.743 

59.34 

20,85 

13.39 

.13094 

6.743 

59.34 

20.85 

13.59 

.13094 

34.59 

59.54 

71.65 

13.59 

.13094 

34.59 

59-3'+ 

71.65 

6,61 

.11330 

34.59 

59.54 

204.85 

6.61 

.11350 

34.59 

59.54 

204.85 

6.61 

.11330 

8,552 

59.34 

206.35 

6.61 

.11330 

8.352 

59.34 

206.35 

1.50 

.00895 

8.352 

209.00 

207.85 

1.50 

.00895 

8.352 

209.00 

207.85 

1.50 

.00895 

94,92 

209.00 

209.95 

1.50 

.00895 

94.92 

209.00 

209.95 

1.50 

.00895 

6.914 

209.00 

211.45 

1.50 

.00895 

6.914 

209.00 

211.45 

.35 

.00895 

6.914 

26.46 

212.95 

.34 

.00895 

6.914 

26.00 

212.95 

.34 

.00895 

13.900 

26.00 

217.95 

.31 

.00895 

10.400 

21.00 

222.95 

.20 

.00895 

8.000 

17.00 

227.95 

.12 

.00895 

6.000 

l4.8o 

252.95 

.05 

.00895 

3.200 

15.60 

256.55 

.047 

.00895 

2.100 

13.10 

236.55 

.047 

.00895 

1.403 

13.10 

237.85 

.045 

.00895 

1.403 

12.92 

237.85 

i.4o 

.01474 

1.403 

51.85 

259.35 

1.393 

.01474 

1.403 

36.89 

259.55 

1.395 

.01474 

4.000 

36.89 

245,25 

1.536 

.01474 

3.570 

56.71 

245.25 

1.336 

.01474 

1.042 

56.71 

246.75 

1.321 

.01474 

1,042 

61.75 

246.75 

1.321 

.01474 

1.042 

61.75 

248,25 

1.516 

.01474 

1.042 

61,75 

248.25 

1,316 

.01474 

3.561 

61.75 

257-85 

1.200 

.01474 

3.361 

61.75 

257.85 

1.200 

.03673 

5.361 

61.75 

258.75 

1.200 

.05675 

3.361 

61.75 

258.75 

1.200 

.03673 

4.910 

33.69 

262.85 

1.200 

.05673 

4.910 

33.69 

262.85 

.385 

.026668 

4.910 

33.69 

424.85 

.385 

.026668 

4.910 

33.69 

424.85 

.385 

.026668 

1.665 

53.69 

426.55 

.385 

.026668 

1,665 

33.69 

426.35 

1.090 

.011718 

1.665 

24.88 

427.85 

1.090 

.011718 

1.665 

24.88 









TABLE I.- PHYSICAL CHARACTERISTICS OF A MULTISTAGE LAUNCH VEHICLE - Concluded 










































































"The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl¬ 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities ay^d the results thereof.^' 

—National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 
of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri¬ 
bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con¬ 
nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 
language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 
and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 
NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availahilify of these publications may be obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 


